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The purpose of this paper is to make an explicit construction of specific self-adjoint 
extensions of the Dirac Hamiltonian in the presence of a <5-sphere interaction of 
finite radius. The exact resolvent kernel of the free Dirac operator is given. This 
specifies related results that have recently appeared in the literature. 



1. Introduction 

Solvable Hamiltonians with singular (5-sphere potentials have been the 
subject of many investigations in the past since the work by Green and 
Moszkowsky.^ These models belong to the small class of models that are 
analytically solvable. As far as we know, all these studies have used the 
von Neumann formalism of self-adjoint extensions of symmetric linear oper- 
ators. The earlier works concern nonrelativistic quantum mechanics based 
on the Schrodinger Hamiltonian with S interactions. A good summary of 
these nonrelativistic models is to be found in Refs. 2 and 3. Comparatively, 
a much smaller number of papers have dealt with the relativistic J-Dirac 
Hamiltonian. To the best of our knowledge, the first rigorous treatments of 
singular Dirac Hamiltonians were provided by Dittrich, Exner and Seba.^'^ 
Next, in the same vein, there have been the papers by Avossevou and 
Hounkonnou,^'^'*'^'^° and more recently those by Shabani and Vyabandi^^ 
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which include quite an uncommon critique of the work in Refs. 6 10. In 
this contribution, we provide self-adjoint extensions of the i5-Dirac Hamil- 
tonian with support on a sphere of finite radius, given specific boundary 
conditions. 

The paper is organized as follows. We present the model and give the 
normalised solutions for the free Dirac radial operator in Sees. 2 and 3. 
This leads us to the construction of the resolvent of the free Dirac radial 
Hamiltonian in Sec. 4. In Sees. 5 and 6 we then construct the resolvent of 
the (5-Dirac Hamiltonian given boundary conditions of the first type. 



2. Model and Boundary Conditions 

Let us consider in the Hilbert space H the formal expression describing the 
interaction model 

H = Hd + V{\x\) , xeW" , (1) 

where is the interaction potential to be specified hereafter. The 

quantity Hd is the free Dirac operator, 

Hd = -ihcaV + Pmc^ , h = /i/27r , (2) 

where the following notations are used, 

-Co) . i-CoX)' 

a = a'' (i- = 1, 2, 3) being the 2x2 Pauli matrices defined by 



^loy' \i J ' vo-i 

so that the Hamiltonian is explicitly given by 

The method adopted for the study of this Hamiltonian is based on 
the von Neuman theory of self-adjoint extensions of linear symmetric 
operators. ^^'^^ The first step is to consider within a Hilbert space H the 
symmetric operator H = Hd, V{H) = {ijj e H^-^{M?) ® , i}){Sn) = 0}, 
where Sr is the closed ball of radius R centered at the origin in R^, and 
iJ^'^ the Sobolev space of indices (1,2). The closure of this operator is 
given by 



H - HD\c^{m?\Sa) - i?£)|jji,2(g3\5^) • (6) 
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The Hilbert space Ti. decomposes as^'^ 

oo J + (l/2) 

W= t^/ [L2((0,cx)),dr)®C2] (g) ,^],,,] , 

j=(l/2) ^=j-(l/2) 

(7) 

where the isomorphism Uji has been introduced to remove the weight factor 
from the measure on the original radial space. The corresponding radial 
self-adjoint operator writes as 

oo i+(l/2) 
j=(l/2) «=j-(l/2) 

= jv' f ■^(''J ) e L2((0, oo), dr) C^; V' e Aaoe((0, oo) \ {R}); 



/ij-^V e i^((0,oo))(g)C2} , (9) 

where /(r) and g(r) determine the radial dependence of the upper and lower 
Dirac bi-spinors. The symbol 1 in Eq. (8) stands for the unit operator. 

Note that the radial operator obtained in Eq. (8) is the same as the one 
given in Eq. (2.9a) of Ref. 4 provided a choice of units such that h=\ = c 
is made. 

The interaction potential associated to this bi-spinor decomposition 
writes as,* 

V{\x\) = G8{\x\-R) = G8{r-R) , G=(^'Jj'^J ,ao,Po&^- (10) 

The boundary conditions we use to characterize the self-adjoint exten- 
sions are obtained as follows. Assume that the limits 'tl>{R±) exist and 
consider the eigenvalue equation 

hji + G5{r - i?)] V = -EV' • (11) 

Then, integrating over {R — e,R + s) taking the limit e ^ 0, and using the 
relation 

fR+e 1 

/ dr6{r-R)i;{r) = -mR+) + ^{R-)] , e^O, (12) 
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we obtain 

Qg + hCTo^ ^{R+) + Qg - f>^T)j ^{R.) =0 , To=(^^'^^ . (13) 

Proposition: The general boundary conditions Eq. (13) define a self- 
adjoint extension of hji iff G+ = G. 

Proof. The proof proceeds in the same manner as that of Proposition 4.1 
in Rcf. 4. 

To solve the problem, henceforth we shall work in a given total an- 
gular momentum sector {j€) of the Dirac spectrum, since the considered 
scalar potential is spherically symmetric. Hence, we shall use the standard 
notation as follows, 

i = j±\ , K = Kii = {-iy-'+'/\j+\) , i = l/2,3/2,5/2,... . (14) 

In what follows, we shall give expressions valid whether for £ = j -|- 1/2 or 
for t = j — 1/2, the upper/lower signs corresponding to these two situations 
in the same order, namely ^ = j ± 1/2. Thus we have 

(15) 

i = j-l: n = -{j + l) = -{t+l) . 

The strategy for relating appropriate choices of boundary conditions 

which define self-adjoint extensions of the Dirac Hamiltonian to the pa- 
rameters fJ'ijiz) in Krcin's formula (see Eq. (46)) in the general setting is 
clear. It suffices to apply this formula to a set of functions {F{r') G{r') and 
determine which boundary conditions are obeyed, as a function of the choice 
for the quantities fiij{z). Here, we shall not pursue the general analysis, 
but only restrict to the two cases which are called^'* "boundary conditions 
of the first type", namely when either one of the two functions /(r) or g{r) 
is continuous at r = R. while the other is discontinuous at r = i? with a 
discontinuity that is proportional to the value of the continuous function 
at that point. Namely, we shall restrict to the following two cases, 

f{R+) = /(i?_) = f{R) , g{R+) - giR.) = af{R) , (16) 

f{R+)-fiR.)^PgiR) , giR+) = g{R.) = g{R) , (17) 

where a = ao/{hc) and f3 = —f3o/{fic) are arbitrary real parameters. As we 
shall see, since such a choice of boundary conditions may indeed be put into 
one-to-one correspondence with a choice for the Krein parameters fJ>ij{z), 
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these conditions define two separate one-parameter self-adjoint extensions 
of the Dirac Hamiltonian, associated to the (5-sphere scalar potential. 

From now on, the boundary conditions defined by Eq. (16) will be called 
"a-type boundary condition of the first type" and those defined by Eq. (17) 
will be called "/3-type boundary condition of the first type" . Note, as it 
should be, that these boundary conditions are particular cases of boundary 
conditions given in Eq. (13). 

The self-adjoint extension of H corresponding to the specific boundary 
conditions in Eqs. (16) and (17) write, respectively, 



(18) 



with 



T^ihje,a,) = {V' e L^((0,oo)) ^ ; fje e AQoe((0,oo)) ; 



gje&ACioc{iO,cx))\{R}y, 
gji{R+) - gjeiR-) = afji(R); 
V,c«V'eL^((0,oo))(g)C2 } , 



(19) 



and 



(20) 



with 




T 



{V e l2((0, (X))) ® ■ gje e ACiociiO, oo)) ; 

fjeGACioc{{0,cx^)\{R}y, 

fje{z; R+) - fji{z; R-) = (3jigji{z; R); 

/ij^,ft-,V'ei^((O,oo))0C2 } . 



(21) 



3. Solutions to the Free Radial Dirac Hamiltonian 



In the present section, we are interested in the determination of the cigcn- 
states of the free radial operator, and in particular their absolute normali- 
sation as required for the construction of the associated free Dirac resolvent 
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or kernel of the corresponding Green's function. Hence, we must solve for 
the eigenvalue equation 



mc 



,2 



f^{-fr + T)\(nr)\_^(f{r) 



The choices of normalisation and signs for these components is standard 
and we required only that 

r dr{r{r;p) g*{T;p)) {^^p^'}.) =S{p-p') . (23) 
Jo \9v^P)/ 

We have to distinguish the two cases of either positive or negative energy 
states, since we know that the free Dirac energy spectrum is that of positive 
and negative energy states such that 

E = ± V(pc)2 + (toc2)2 , p>0 . (24) 

Then, introducing the following parametrisations, 



E{p) = ±u;{p) , w{p) = V(pc)2 + (mc2)2 , p=^r, (25) 

and labelling all distinct eigenstates by the value for p > and the sign of 
the energy, one finds for the normalised solutions to the free radial Dirac 
Hamiltonian of positive energy, for £ = j ± 1/2, 

^^^''^'^ = V2!^ip)iX)-mc^) ''^'^^ ' ^''^ 

The negative energy solutions are given by, with £ = j ±1/2, 

^-^'''^ = V2^(p)(r(p) + mc2) (ir)'''^^-^/^ (¥) ' 

, , —pc (pr\^l'^ (VT\ 

Here, Jp(-) arc Besscl functions of order [i}^ Finally, let us note that one 
may readily check that the configurations of opposite energy have a van- 
ishing inner product, given the negative sign in the function g_{r;p) for 
the negative energy solutions relative to that of the positive energy ones. 
Hence, we have constriicted a complete basis of eigenstates for the free ra- 
dial Dirac Hamiltonian, in terms of which it is thus possible to construct 
the explicit expressions for the associated resolvent. 
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4. The Resolvent of the Free Dirac Radial Hamiltonian 

In an abstract and formal manner, let us consider a self-adjoint operator A 
with spectrum* 

^IV-n >= A„|V'n > , (30) 

where the eigenstates are orthonormalised, 

< i>n\tlJm >= Sn,m ■ (31) 

Consequently, given a complex parameter z, the associated resolvent 
[A — z] is simply given by the representation 

[A-z]-'=J2\^n>Y^<'l'n\, (32) 

where z does not belong to the >l-spectrum. We shall simply follow exactly 
the same construction for the (r, r') kernel of the free Dirac radial Hamil- 
tonian, since we have constructed its spectrum in the previous section. In 
our case, this resolvent is also a 2x2 matrix, each of which elements are 
obtained in a likewise manner. Hence, we have for the resolvent kernel of 
the free Dirac radial Hamiltonian 



Go(r,r';.)=(-"^l'': • (33) 



Gn{r,r';z) Gi2(r,r';z) 
G2i{r,r';z) G22{r,r';z) 

Introducing the quantity k = k{z) defined by 

k{z)c = _ (^c2)2 , (34) 

with the usual branch cut along the negative real axis for the square-root 
function, we have for the Gu element 

i-KZ + rru? ( k[z)c\ f Hz)c, 



Gniry;z) = -^^Vr^J,,,^2[^r) i^-/, ^-^r' ) , (36) 

Eqs. (35) and (36) being valid for r > r' and r < r' , respectively. Note that 
as a consequence, the quantity Gii(r,r' ; z) is also well defined for r — r' , 
as it should. 



''We assume here that the spectrum is not degenerate, since this is the situation encoun- 
tered presently. The case of a degenerate spectrum is a straightforward generalisation. 
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The second diagonal element is given by 

, inz-mc^ , — ^ (i) f k{z)c \ ( k{z)c \ 

(37) 

, mz-mc^ J—; (k{z)c \ (i) ( k{z)c ^ 

(38) 

Eqs. (37) and (38) being valid for r > r' and r <r', respectively. 

For the off-diagonal element Gi2{r,r'-, z), we obtain the following two 
expressions valid for r > r' and r < r' , respectively, 

= ?;g|v^<.(^0 w.(^/) . m 

e„(,,;.,.^«|v^7.„,.(M|£r) <U.(ag;.) . ,40, 

Finally for G2i{r, r'; z) element, we obtain 

Eqs. (41) and (42) being valid for r > r' and r < r' , respectively. 
5. The Deficient Index Subspaces 

In order to use Krein's formula to express the Dirac radial resolvent for a 
given choice of boundary conditions which define a self-adjoint extension 
of the Dirac Hamiltonian related to a 5-sphere scalar potential, it is first 
necessary to establish a basis for the deficiency subspaces. This amounts to 
identifying the cigcnsolutions to the previous eigenvalue problem in which 
the energy eigenvalue E is now replaced by a complex parameter z which 
does not belong to the spectrum of the free and interacting Hamiltonians. 
The only restriction is that the obtained solutions must be normalisable for 
the norm dr. 

Consequently, for a given complex value of z, let us introduce again the 
quantity 



k{z)c = _ („c2)2 , (43) 
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where the square-root function is defined as previously. It is then straight- 
forward to show that there exist two unique Unearly independent normal- 
isable solutions, for ImA;(2;) > 0. The first solution is given by, 



(44) 



.if r>i?: /i(r; 2) =0, 
g\{r\z) = . 

The second solution is given by 

.ifr<i?:/2(r;2)=0, 
giir-.z) = ; 



. u r > K . f2[r, z) - ^^———^^ y rj u^^^^^ \ "^J ' 

If lmk{z) < 0, the corresponding deficiency subspace is spanned by 
the complex conjugates of the above expressions, except for the fact that 
k{z) is then still given by the above definition in terms of z rather than in 
terms of z* = z. In other words, for Imfc(2:) < 0, the deficiency subspace is 
spanned by the same solutions {fi{r; z) gi{r; z)) as above and by functions 
(/2(r;2;) 32(^5-^)) which are given by the same expressions as those for 
if2{r; z) g2{r; z)) with the function -ff£+\/2=Fi ^'^P^^'Ced by -ff]+i/2^i- 

Since Krein's formula reads^^ 



G(r,r';z) = Go(r,r';z) + ^ /..,(z) (//(r'l^) 



0) , 



(46) 

where G(r, r'; z) stands for the resolvent kernel of the interacting Hamilto- 
nian, it is clear that the specific values for the coefficients fjbij [z) are defined 
up to a normalisation which is correlated to the chosen normalisation for 
the solutions (/j gi) constructed above. As mentioned already, the lat- 
ter normalisation is simply chosen by analogy with that of the free energy 
eigenstates, and it might thus turn out that there could be a choice of 
normalisation rendering the expressions for the coefficients fJ-ij{z) simpler. 
This is indeed the case, but in a manner which is dependent on the specific 
choice of boundary conditions which is associated to a specific self-adjoint 
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extension of the Dirac Hamiltonian. It is for this reason that we have kept 
the above normahsation, leaving it as an exercise to change the normal- 
isation to make it simpler once a specific choice of self-adjoint boundary 
conditions is considered. This is rather straightforward on the basis of the 
results to be presented in the next section. 

Note also that since the deficiency indices of the Hamiltonian arc (2, 2), 
the set of its self-adjoint extensions is characterized by four independent 
real parameters, which arc thus in correspondence with the four param- 
eters Hij{z). In other words, any choice of boundary conditions for the 
functions /(r) and g{r) obeying the interacting Hamiltonian energy eigen- 
value problem which may be put into unique correspondence with the pa- 
rameters ^iij{z) in Krein's formula determines a self-adjoint extension of 
the Dirac Hamiltonian. The next section displays two such choices, each 
characterized by a single real parameter, leaving it aside how to construct 
the general case and how to identify its relationship to 5-sphere scalar and 
vector interactions. 



6. The Resolvent Equation 

6.1. The OL-Type Boundary Condition of the First Type 
Requiring /(r) to be continuous &t r = R implies that we must have 

/ \ I \tt(i) /k{z)R. , , . k{z)R 

iiii[z) ^ Hi{z)Hl^^,^{ ) , H2i[z) ^ iii{z)Je+i/2[^—) , (47) 



Ml2(^)=M2(^)if,%2(^) , M22(^) = M2(^)Jm/2(^), (48) 

where i-ii{z) and ii.2{z) arc real parameters still to be determined. After 
some algebra, one finds that the boundary condition involving the difference 
g{R+) — g{R-) implies the following relations for the quantities fJbi{z) and 
IJi2{z) in terms of a, for (. = j ± 1/2, 

. . iiT R(z + mc^) in 2hz 1 



hc 2 he k{z)c 
-(1) 

^-1-1/2 



„(1) ,kj^. 
"-"^-1-1 /2V ft ) 



" 2 he ^i+i/2\ h j'^l+^/'^K h ) 



(49) 
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IJ,2iz) = 



in R{z + mc^) in 2hz 1 



he 



2 he k{z)c 



a J, 



1+1/2 



(50) 



he 



Note that the quantity that multiphes a in the denominator of these 
two expressions is simply the value for Gii{R,R;z) for the free Dirac re- 
solvent, in exactly the same manner as for the nonrclativistic Schrodinger 
problem. Furthermore, it should be clear that most of the overall nor- 
malisation factors multiplying the last fraction could be absorbed directly 
through a renormalisation of the functions fi{r] z) and gi{r\ z) spanning the 
deficiency subspacc. Wc shall leave this point open, since this is rather triv- 
ial, even though it simplifies somewhat the expression for the coefficients 
IJL\{z) and lJL2{z). Finally, notice that the final expressions for all four coef- 
ficients fJ-ij{z) generalise similar expressions relevant to the nonrclativistic 
Schrodinger problem. In the latter case, all the factors involving the Bessel 
functions evaluated at k{z)R/h may be absorbed into the normalisation of 
the vectors spanning the deficiency subspace, leaving over a rather simple 
expression for the coefficient /.j,{z) in that case. For the Dirac Hamiltonian, 
this is no longer the case in a simple way, but still the general structure is 
maintained for the solutions of the coefficients lJ-ij{z). 

Finally the resolvent of the extended Dirac Hamiltonian for the a-type 
boundary conditions is given by 



{hi. 



{hi,o - z)-^ + A{z)l2i^{r; z) {r{r;z), •) , (51) 



with the following notations, 

in R{z + mc^) in 2hz 1 



A{z) = 



he 



2 he k{z)e 
a 



-rl Rjz+mc^) tt(1) ( k{z)R ,^j j k{z)R ^ 



(52) 



1/^1 (r;z) 
fp2{r;z) 



(53) 



ipi{r;z) = a{z) x ( ^^^r 



J, 



e+i/2 



h 



+ Je.+i/2{ 



><z)R (1) Mzy 



h 



)^e+i/2i 



-) 



(54) 
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ip2{r;z) = b{z) x 



Je+i/2Ti{^—)^e+i/2y^—)+ 



Mz)R ^i) k{z)r 



(55) 



a{z) = ^Jk}^= , b{z) = ^Jk}^=. (56) 

^/2hz{z — mc^) ^2hz{z + 'mc^) 

6.2. The (3 -Type Boundary Condition of the First Type 

A similar analysis is possible when the continuity requirement is imposed 
for g{r), 

g{R+) = g{R.) = g{R) . (57) 

It then turns out that the /3-type boundary conditions of the first type, 
Eq. (17), are then possible, and provide thus a one-parameter self-adjoint 
extension of the Dirac Hamiltonian. In this case, we must have 

, . , ,.,(1) ,k(z)R. , , , . ^ ,k{z)R. 

Mii(^) =Mi(^)i^£+Wi(— ^) ' ^^2i[z) = ni{z)Ji+i/2^i{-^—) , (58) 

Ml2(-^)=M2(^)<\/2^l(^) , ^i22{z)=^,2{z)Je+„2^,{^^), (59) 

while the coefficients and (i2{z) are given as follows in terms of the 

real parameter /3 defining the boundary conditions Eq. (17), for £ = 1/2, 

i-K R{z — mc?) in 2hz 1 



fj,i{z) = 



he 2 he k{z)e 



(60) 



H2{z) = 



in R{z — mc?) in 2hz 1 



he 2 he k{z)e 

h 



PJi+l/2^lC-^' 



Qiir R(z-mc^) ujl) / k{z)R -^ j / k(z)R. 

~ hi ^^-Hl/2Tl*-~fi~''"'^+l/2TU— 



(61) 



The same remarks as those made for the a-type boundary conditions 
are of application in this case, of course, with in particular the value for 
G22{R,R;z) explicitly multiplying the parameter (3 in the denominator of 
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the last fraction factor. Finally, the resolvent in a compact form, for the 
/3-type boundary conditions, is given by 

{h,j,, - z)-' = {hi,o - z)-' + A{z)l2Hr; z) {r{r;z), •) , (62) 

with the notations 

iJi{r]z)\ 
i^2{r;z) I 



il){r; z) 



(63) 



V'i(r;z) = a{z) x 



1/2 r 



'Je+i/2(^—)^e+i/2Ti^ )+ 



k{z)R (^i) ,k{zy. 

+ -'m/2=Fl(— ^)-«^+l/2(— ^) 



(64) 



'fc(z) ^ 
ij^2{r;z) = b{z) x ( -—r 



J f^^^^'' m^^) Mz)R . 



k{z)R >i) , Hz)r , 

+ Ji+l/2:fl{^—)^e+i/2^:l{^—) 



(65) 



~, . iir R{z — mc^) in 2hz 1 

^^'^ Y T he k{z)c ^ 



(66) 



7. Conclusion 

We have thus constructed the resolvent kernel G{r,r';z) for the one- 
parameter self-adjoint extensions of the Dirac Hamiltonian associated to 
the ^-sphere scalar potential and characterized by the a- and /?-type bound- 
ary conditions of the first type defined in Eqs. (16) and (17). More general 
cases could be considered on the basis of the general expressions obtained 
above for the values of f{R±) and g{R±), in order to identify all possible 
self-adjoint extensions parametrized by four independent real parameters. 

Associated to the boundary conditions Eqs. (16) and (17), these results 
should now enable the analysis of the scattering, spectral and resonance 
properties of the constructed self-adjoint extensions of the Dirac operator.^"' 

Finally, let us say that it is also possible to construct a subclass of self- 
adjoint extensions characterized by two parameters. Shabani and Vyabandi 
attempted such a construction in their work^^ but it is not clear to us how 
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the resolvent kernel of the free Dirac operator whieh they obtained and on 
which the remainder of their analysis rests in a crucial way, compares to its 
explicit contruction as detailed in the present contribution. 
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